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ABSTRACT. New exact solutions for unsteady magneto-hydro dy- 
namic (MHD) flows of a generalized second-grade fluid due to a 
impulsive motion of the plate have been derived. The generalized 
second-grade fluid saturates the porous space. Fractional deriva- 
tive is used in the governing equation. The analytical expressions 
for velocity and shear stress fields have been obtained by using 
Laplace transform technique for the fractional calculus. The ob- 
tained solutions are expressed in series form in terms of Fox H- 
functions. 


1. INTRODUCTION 


During the past several decades, the Navier-Stokes equations 
were intensively studied in the literature. However, these equations 
are unable to predict the behaviors of many real fluids encountered in 
numerous industrial applications such as processing of polymers, phar- 
maceuticals, personal care products, food products, inks, and so forth. 
Therefore, it has been now recognized in industrial and technologi- 
cal applications that non-Newtonian fluids are more appropriate than 
Newtonian fluids. Non-Newtonian fluids are a broad class of fluids in 
which the relation connecting the shear stress and shear rate is non 
linear and hence there is no universal constitutive model available that 
can alone predict the behavior of all non-Newtonian fluids. Hence, it 
is necessary to study the behavior of non-Newtonian fluids in order 
to obtain a thorough understanding and improve their utilization in 
various manufacturing processes. Rivlin and Ericksen [1] introduced 
a subclass of non-Newtonian fluids known as second-grade fluid for 
which a possibility exist to obtain the exact solution. Exact solutions 
of second-grade fluid for start-up flows have been investigated by Ban- 
delli [2] using integral transform technique. Tan [3] discussed the flow 
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of suddenly moved flat plate in a generalized second-grade. Exact solu- 
tions of a generalized second grade fluid corresponding to the oscillatory 
flow between two cylinders have been achieved by Mahmood et al. [4]. 
Tripathy [5] discussed peristaltic motion of a generalized second grade 
fluid passing through a cylindrical tube. Tan [6] obtained solutions for 
unsteady motions between two parallel plates of the generalized second 
grade fluid. 

In the last few decades the study of fluid motions through porous medi- 
um have received much attention due to its importance not only to the 
field of academic but also to the industry. Such motions have many 
applications in many industrial and biological processes such as food 
industry, irrigation problems, oil exploitation, motion of blood in the 
cardiovascular system [7], chemistry and bio-engineering, soap and cel- 
lulose solutions and in biophysical sciences where the human lungs are 
considered as a porous layer. etc. Unsteady MHD flows of viscoelastic 
fluids passing through porous space are of considerable interest. In the 
last few years alot of work has been done on MHD flow, see [8-12] and 
reference therein. 

Recently, the fractional derivative [13] approach is proving to be an im- 
portant tool for considering the behaviors of such types of fluids. Many 
researchers investigated different problems using fractional derivative 
technique regarding such fluids. In their works, the integer order time 
derivatives in the constitutive models for generalized second-grade flu- 
ids are replaced by the Riemann-Liouville fractional derivatives. Alot 
of work has been done on fractional derivatives during the last few 
years. Here we mention only those contributions which regards with 
the viscoelastic type fluids [13-20] and the references therein. 
According to the authors informations upto yet no study has been done 
on the MHD flow of generalized second-grade fluid induced by impul- 
sive motion of the plate flowing through a porous space. Hence, our 
main objective in this note is to make a contribution in this regard. 
We take an incompressible MHD flow passing through porous space of 
a generalized second-grade fluid. Laplace transform method has been 
used for the fractional calculus to obtained exact solutions for the pro- 
files of velocity field and the corresponding shear stress. The obtained 
solutions satisfies all the imposed initial and boundary conditions are 
expressed in terms of Fox-H function. 
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2. GOVERNING EQUATIONS 


The equation of continuity and momentum of MHD flow passing 
through porous space is given by:(Tan and Masuoka [6]) 


(2.1) V.V-0 e) = divT — of? V +R, 

where V=(u,v,w) represents velocity vector, electrical conductivity and 
density of the fluid are represented by c and p respectively, Bo is the 
magnitude of a uniform magnetic field, material time derivative is de- 
noted by d/dt, Cauchy stress tensor is represented by T, and R is the 
Darcy's resistance of the porous space. 

For an incompressible and unsteady generalized second-grade fluid the 
cauchy stress tensor T is given as [7]: 


(2.2) T —S-pL S-—y,uW,;-o1Ws- o03W3, 


where S and pl represents the extra stress tensor and the indeterminate 
spherical stress, the dynamic viscosity is denoted by jz, normal stress 
moduli are represented by a, and a2 and the kinematic tensors are W, 
and W^» defined as 


23 W, =L+L’, W = D? - WIL- L^W, 
t 


where L is the velocity gradient and D? represents the operator for 
fractional differentiation whose order is @ and is based on the Riemann- 
Liouville definition [13], 





1 d (^ g(t) 

24 D? = dt, 0<b<1 

G4) — Dil = sya | ope os 
where Gamma function is denoted by I'(-). Model for ordinary second- 
grade fluid can be obtained by putting 6 — 1. For the compatibility 
of this model with thermodynamics it is required that the material 
moduli should obey the following conditions 

(2.5) Q4 T @ = 0, 0 € o4 and u > 0. 

For the second-grade fluid the Darcy's resistance satisfies the following 
equation: 

@ o 

2.6 R=-—-+-(1 —)V 

(2.6) Elta) 
where k > 0 and 9(0 < ¢ < 1) are the permeability and the porosity of 
the porous medium. For the following problem we consider the velocity 
field and an extra stress of the form 


(2.7) V = (u(y, t),0,0), S = S(y,t). 
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where u is the velocity taken in the x-direction. Substituting Eq.(7) 
into Eq.(2) and taking into account the initial condition 


(2.8) S(y,0) — 0, y > 0, 
the fluid being at rest up to the time t = 0, we get 


(2.9) Sey = (u + o1 Dj )Oyuly, t), 





where Syy Sue Sus Syz 0, and Sz, = Sys- The balance of 
linear momentum in the absence of body forces and pressure gradient 
is given as: 





Lp 0 
pus + oix) T pojuly, t), 


(2.10) Oy Sey — o Bau — 
By putting S,, from Eq. (9) into (10), we find the governing equation 
under the form 

(2.11) 


o 
pðiuly,t) = (u + o D puty, t) — o Buly, t) — (iz). 


3. STATEMENT OF THE PROBLEM 


We take an unsteady incompressible flow of homogenous and electri- 
cally conducting second-grade fluid bounded by a rigid plate at y = 0. 
The plate is taken normal to y-axis and the fluid saturates the porous 
medium y > 0. The electrically conducting fluid is stressed by a unifor- 
m magnetic field Bọ parallel to the y axis, while the induced magnetic 
field is neglected by choosing a small magnetic Reynolds number. Ini- 
tially, both the plate and the fluid are at rest, and after time t=0, it is 
suddenly set into motion by translating the flate plate in its plane, with 
a constant velocity A. The initial and boundary conditions of velocity 
field are: 


u(y, 0) = 0; y > 0, 
(3.1) HODES: 


u(y,t), Oyu(y,t) +0 as y — oo and t > 0. 
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4. CALCULATION OF VELOCITY FIELD 


Employing the non-dimensional quantities 


u * yU * tU? * aU? * A 
T5 uy Bu t d Q 7 ,A E 
(4.1) 





_ sS — RU? 2. cvBj 
esp K= grm M = 202 > 





T 


The dimensionless mark * is omitted here for simplicity. Thus, the 
governing equations of dimensionless motion becomes 


(4.2) Opu(y,t) = (1+ aD) )O2uly, t)— La + oS uly, t) — M?u(y, t), 


(4.3) r(y,t) = (1+ aD?)O,u(y, t) 


with the given conditions as 


u(y,0) = 0; y > 0, 


i) HOD SATS 


uly, t), Qyu(y, t) +0 as y — oo, and t > 0. 


First we will apply the Laplace transform to eq (14) and using the 
Laplace transform formula for sequential fractional derivatives [21] 


(4.5) u(y, t) a u(y, te “dt, s > 0, 


Taking into the account the corresponding initial and boundary condi- 
tions (16), we get the following differential equation 


z 1+aq q+ M? 
4. 9? - | 
uo gawd- (at 








» A 
u(0,g)-— —; t2 0, 
q 
(4.7) u(y,q), Oyu(y, q) —> 0 as y > oo, and q > 0. 


The solution of Eq.(18) satisfying the boundary conditions (19) is of 
the following form: 





11 A 1 l | | 2 
(4.8) u(y,q) — EET Faq) + Kiq + M » 
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To get the analytical solution for velocity field and to avoid difficult 
calculations of contour integrals and residues, we will apply the discrete 
inverse Laplace transform method [21], but first we have to expressed 
Eq. (20) in series form as 





(4.9) 
eitfitgithi tT +S cyhi tr s M?n38 


E SS ER EUR LUN Pes sA A(-1) 
23600 NDS ey filgilhaylrlsUcei/2-ig- fi-i-Br-s 


e1=0 fi-0 g1—0 h4—0 r=0 s=0 
I(fh —e/2)F(g1 — fa) (a + f(r + e/2)(s — ei/2) 
T'(e1/2)T (7e1/2) (e3/2)P CR) P (A) 
Now apply the discrete inverse Laplace transform to Eq. (22), we get 


























(4.10) 
yertfitnthitrtsy fi-hi—Br Sqhitrt+s 


Hon UR SUE %2 Ser es A(-1 
u(y, t) = 2: p» 2s X 2 > ey filglhilr!s T (— fi) Ke12-f 


e1=0 f1=0 91-0 hi=0 r=0 s=0 
T(f; — e/2)T (gi — fa (ha + fa) P(r + ei/2)T (s — €1/2) M? ye 
L0e3/2)EC—e172)V Ces 2) 3) (fı — hy — BP) 


To get Eq. (23) in a more compact form we use Fox H-function |13], 

















(4.11) 
œo oo 00 eit fiegirts yp2g1,01t—fi-Br-stlgrts 


B UR um MEC. (—1) 
up poo — eifilgirisKei2- 
(1— fı +e1/2,0), (1 — gı + fı1,0), o 1), (1 — s + e1/2,0), 
a (1—r —e1/2,0). 
5.7 t (1 m &1/2,0), (1 = Tis Oh (1 + fi, 0), (0, 1), (1 ns e3/2, 0), 
(1— €1/2,0), (fi + Br + 5, —1). 


To obtain Eq. (24), the following Fox H-function property has been 
used, 


l,s 
H | 


s,t+1 

















(1 — ay, Aj), ..., (1 — as, As) | B T D(a4 + Ar)... (as + A,r) , 


(1,0), (1 — bi, Bi), o (1755 Be) | ~ Ser + Bir). TG F Br)” 


5. CALCULATION OF SHEAR STRESS 


To get the shear stress first we apply Laplace transform on Eq. (15), 
we get 


(5.1) T(y,q) = (1 + aq”) dyui(y, q), 
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Substituting u(y,q) from eq. (20), we get 
A(1 P 
(5.2) T(y, t) = -AL ta ag - V By) VB. 


where 








To get a more compact form of 7(y, q), we write eq. (26) in series form 
as 


(5.3) 


ja +fitgithitd+r+s+l 


=> x 3 S ON LL OA SGS 


e1=0 f1=0 g1=0 hi=0 r=0 s 
qhiththtmitr+sp( fi — eii/2)P (gi — 4 (^i + fi) (tr + e&/2)T (s — e1/2) 
T'(e/2)T(—e:/2) (3/2) CF) ( - A) K275i 
Mn+ AyD (i — 1/2) (ji — 4)T (ki + 1) (l5 — 1/2) (m; — 1/2) 
l'(1/2)T (1/2)T (1/2) (6) (tag 75 7h Bh m Brent 














x 











where 





ok OO œ © © oo 
^ —^ VA NA NA v 
a 0 j1 0 ky Ol, 0m =0 

Ġ =i +j +ki+th +m, 











Taking the inverse Laplace of eq.(27), we get 


(5.4) 


jed fitgithitatr+stl 


257 » » Y 25 EX DEDI 


e1=0 fi =0 g1—0 hi=0 r=0 s 
alra Pareres p fy — e,/2)0 (gi — fl (a + fil (r + e1/2)T (3 — €1/2) 
T(ex/2)F(—e1/2)L (1/2) T (FT (f) PR en 
M?s?5,6T(i, — 1/2) T (j; — 4)T (s + 4)E (I — 1/2) (m4 — 1/2) 

















x 








PZPP UDT PA — hi — 1 — k — BL, — my — Br - 8 +1) 


x t fi-hi—i1—-ki—Bly-—m1—Br—s 
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Finally, using Fox H-function to get the stress field as, 


(5.5) 


jen +fitgitG trts+iyer 


ay 3 y 25 x WEG HP 29 25i 


e1=0 f1-0 gi=0 r=0 s=0 
p ij—ki—8li—mi-—8r—s 

















x 








a kı—lı—-mı-r s Ke/2 fi—à 1/2 
(—i, 4-3/2,0),(1 — jı + 71,0), (1 — kı — 45,0), (1 — h + 1/2,0), 
sae (1 — fi,1),(1 — s + e,/2,0), (1 — r — e1/2,0), (1 — fı + e3/2,0), 
"| (—m + 3/2,0), (1 — 91 + fi, 9). 
E (1/2, 0) (T — i1, 0); (1 + i1, 0), (1 — 61/2, 0); (1 — fı, 0), (0, 1), 
(1+ fi,0), (1 — e1/2, 0), (1/2, 0), (1/2, 0), (1 + e1/2, 0) 
(fi tis ky 8h +m + Br +s, —1). 
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